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ABS TRACT 

A  three  dimensional  finite  element  code  is  written  for 
truss  analysis  and  design.   Trusses  may  be  designed  for  mini- 
mum weight  subject  to  constraints  on:   member  stresses,  Euler 
buckling,  joint  displacements  and  system  natural  frequency. 
The  optimum  configuration  may  be  found  in  addition  to  optimi- 
zation with  respect  to  member  sizes. 

The  finite  element  code  may  be  used  as  a  stand  alone 
analysis  tool  or  may  be  coupled  to  an  optimizer  of  the  user's 
choice.   The  finite  element  displacement  method  of  analysis 
is  used  for  static  analysis  and  eigenvalues  are  calculated 
using  the  subspace  iteration  technique. 

A  major  portion  of  this  work  is  the  software  (user  guide) 
which  is  presented  in  detail  with  examples  and  results .   Ex- 
planation of  how  this  code  is  coupled  to  an  optimizer  is 
given. 

Design  variables  are  member  area  sizes,  joint  coordinates, 
or  both.   Examples  are  presented  to  demonstrate  the  method. 
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I.   INTRODUCTION 

Design  optimization  requires  the  minimization  or  maximi- 
zation of  some  parameter.   Optimization  of  structures  has  had 
continuing  changes  since  its  development  in  the  early  1960's 
with  an  active  area  of  research  being  elastic  truss  structures 
The  main  goal  is  to  design  structural  systems  that  efficiently 
perform  specified  purposes.   For  any  design  to  be  acceptable, 
it  must  satisfy  a  variety  of  physical,  aesthetic  and  economic 
constraints . 

Since  most  physical  problems  can  be  modeled  by  some  mathe- 
matical formulation  a  computer  program  can  be  written  to  per- 
form the  necessary  calculations . 

The  purpose  of  this  research  was  to  develop  a  FINITE 
ELEMENT  code  that  could  be  easily  coupled  to  an  optimizer, 
thus  (1)  allowing  testing  and  comparison  of  various  optimizers 
and  (2)  provide  a  useful  design  program  in  its  own  right. 

The  design  problem  considered  in  this  study  is  the  opti- 
mization of  three-dimensional  indeterminate  trusses ,  for 
multiple  static  load  conditions.   The  objective  is  to  mini- 
mize the  weight  of  the  structure  where  the  design  variables 
are  member  sizes  and  joint  coordinates.   Constraints  include 
stress,  displacement,  Euler  buckling,  and  the  natural  fre- 
quency of  the  system. 
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This  document  describes  the  use  and  capabilities  of  the 
finite  element  computer  code  to  be  coupled  to  an  optimizer. 
The  user's  manual  presented  in  Chapter  V  contains  a  simple 
design  example  in  which  the  program  is  coupled  to  the  CONMIN 
optimization  code  [Ref.  1]. 

Additionally,  guidelines  for  coupling  the  code  to  an 
optimizer  of  the  user's  choice  are  presented. 

Several  examples  demonstrating  the  program  under  a  variety 
of  conditions  are  presented.   Conclusions  and  recommendations 
for  future  work  are  given. 
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II.   OPTIMIZATION 

A.   INTRODUCTION 

The  main  goal  of  structural  engineering  optimization  is 
to  design  structural  systems  that  efficiently  perform  speci- 
fied purposes.   Selection  of  a  specific  algorithm  must  include 
the  following  considerations:   1)  the  structure  should  be 
analyzed  as  few  times  as  possible,  2)  the  algorithm  should 
minimize  specific  gradient  information  required,  and  3)  it 
should  provide  reasonable  assurance  that  an  optimum  design 
will  be  reached. 

The  next  few  sections  discuss  the  general  formulation  of 
the  algorithm  with  respect  to  the  above  requirements . 

3.   FORMULATION 

Minimum  weight  design  of  trusses  is  presented  in  the 
general  form  of  a  mathematical  programming  problem  as  follows: 

Minimize  F(X)  (Eq.  1) 

Subject  to: 

g-  (x)  <_  0  j=l,m        (Eq.  2) 

x-  <  x.  <  x.  i=l,n        (Eq.  3) 

l—i—i  '  H 

where  F(X)  is  the  objective  function,  in  this  case,  weight  of 
the  structure  be  minimized.   X  is  the  vector  of  independent 
design  variables  and  here  contains  the  member  cross-sectional 
areas  as  well  as  the  coordinates  of  the  joints  .   The  inequality 
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constraints,  g.Cx)  <  0,j=l,m,  must  be  satisfied  for  the  design 
to  be  accepted  as  feasible .   These  include  limits  on  stress , 

Euler  buckling,  joint  displacements  and  the  first  fundamental 

o       u 
frequency  of  the  structure.   Side  constraints  x.  and  x.  are 

H     J  11 

lower  and  upper  bounds  on  the  design  variable.   These  may  be 
treated  as  general  inequality  constraints . 


C.   PARAMETERS 

The  design  variables,  constraints,  and  objective  function 
considered  in  the  design  process  are  discussed  here. 

1 .  Design  Variables 

The  design  variables  the  user  has  available  are  member 
cross-sectional  area,  reciprocal  of  member  cross-sectional 
area,  joint  coordinates,  or  both  member  areas  (or  their  reci- 
procals) and  coordinates.   The  cross-sectional  areas  are  A,  , 
k=l,NE,  where  NE  is  the  number  of  elements  in  the  truss 
structure.   The  joint  coordinates  are  X..,  i=l,2,3  j=l,  N J , 
where  i  is  the  coordinate  axis,  j  is  the  joint  number  and  NJ 
is  the  total  number  of  joints.   Design  variable  linking  is 
allowed  in  both  member  sizing  and  coordinate  design  variables. 

2 .  Objective  Function 

The  objective  function  considered  here  is  weight. 

NE    AL 
Weight  (W)  =  Z      p.   i  (Eq.  4) 

i=l   x 

where  p  is  the  material  density,  A  is  the  cross-sectional  area 
and  L-  is  the  length  of  the  member.   The  truss  may  be  made  up 
of  members  of  differing  materials . 
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3.   Constraints 


Stress 


a  .  . 

-il    -l   <   0 


i=l,NE 
j=l,NLC 


(Eq.    5) 


a .  . 

il-l<0 
a 


i  =  1 , ME 
j=l,NLC 


(Eq.    6) 


where  a  is  the  stress  in  member  i  under  load  condition  j .   NE 
is  the  number  of  elements  and  NLC  is  the  number  of  loading 
conditions . 

a        is  the  lower  bound  on  the  stress  (maximum  compressive 
stress)  and  a    is  the  upper  bound  on  the  stress.   The  upper 
and  lower  bounds  on  stress  may  be  different  for  each  member, 
but  are  taken  to  be  the  same  for  every  loading  condition. 

t .   Euler  Buckling 

The  stress  at  which  Euler  buckling  occurs  is  given  by 


a  . 

i 


-K-A.E. 
ill 


i=l,NE 


(Eq.  7) 


where  the  subscript  i  corresponds  to  the  member  number,  E. 
is  Young's  modulus  and  K.  is  a  constant  depending  on  the 
cross-sectional  geometry  of  the  member. 
5 .   Displacement 

Displacement  limits  are  imposed  at  prescribed  joints 
to  create  a  constraint  equation  as  follows: 
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u.  . 


ij^  -1.0  (Eq.  3) 


u.  ., 


u 


ii£  -10  (Eq.  9) 


u .  ., 


where  u.  .."  and  u.  .,  are  lower  and  upper  bounds  on  the  displace' 
ljk       ljk  ^r  v 

ment  at  joint  i  in  the  coordinate  direction  j  under  loading 
condition  k. 

6 .   Frequency 

The  first  fundamental  frequency  of  the  structure  is 
required  to  exceed  the  specified  lower  bound,  so  that 

1  -  X/X    <0  (Eq.  10) 

Reference  2  is  an  excellent  source  for  the  basic 
structural  design  formulation. 
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III.   FINITE  ELEMENT  METHOD 

A.  INTRODUCTION 

Several  features  are  desirable  when  the  finite  element 
methods  of  analysis  are  used  in  design  optimization.   First 
the  number  of  analyses  for  the  structure  should  be  kept  to  a 
minimum.   Second,  the  amount  of  gradient  information  required 
during  the  design  process  should  be  reduced  to  shorten  run 
times  and  computer  storage  requirements.   Third,  the  user 
should  be  able  to  specify  only  that  gradient  information 
desired . 

B.  ANALYSIS 

Initial  problem  formulation  includes  member  sizing,  mate- 
rial properties  (which  may  be  different  for  each  member) ,  a 
set  or  sets  of  external  loads,  and  specified  support  condition. 

The  analysis  for  the  stresses  and  deflections  must  satisfy 
the  conditions  of  equilibrium  of  forces  at  the  nodes  and  com- 
patibility of  deformation.   In  this  analysis  the  weight  of  the 
individual  members  are  not  part  of  the  specified  load  conditions 

Additionally  for  this  analysis  the  following  assumptions 
are  made.  Trusses  will  be  treated  as  discrete  elements,  and 
each  element  will  be  treated  as  pin-connected  with  loads  and 
reactions  supported  at  the  joints. 


16 


The  general  method  of  solution  is  as  follows .   The  Dis- 
placement (Stiffness)   method  reference  2  considers  the  joint 
displacement  components  as  the  unknowns ,  and  written  in 
matrix  notation  for  the  general  case  is 

Ku  =  P  (Eq.  1) 

where  K  is  the  global  stiffness  matrix,  P  is  the  vector  or 
vectors  of  applied  loads,  and  u  is  the  vector  of  vectors  of 
displacements . 

Once  displacements  at  every  node  are  known,  the  internal 
forces  and  stresses  are  calculated  by  applying  the  appropriate 
force-deflection  relationships. 

When  the  system's  natural  frequency  constraints  are  con- 
sidered the  design  process  requires  the  solution  of  an  eigen- 
problem.   This  solution  will  determine  the  natural  frequencies 
and  normal  modes  of  the  structure.   For  linear  elastic  struc- 
tures, the  finite  element  approach  leads  to  the  following 
equation  of  motion  with  free  vibration  conditions, 

Mii  +  Ku  =  0  (Eq.  2) 

where  M  is  the  global  mass  matrix,  and  u  is  the  linear  accel- 
eration vector.   This  leads  to  an  eigenvalue  problem  of  the  form 
1K-X  ld=0,  and  is  solved  here  by  the  subspace  iteration  method 
reference  2 . 
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C.   GRADIENTS  OF  CONSTRAINTS 

Gradient  computation  of  necessary  functions  in  a  design 
optimization  process  arise  from  the  need  for  derivative  infor- 
mation for  efficient  mathematical  programming.   For  structural 
optimization  the  gradient  of  displacement  with  respect  to  the 
design  variables,  9u/3x,  is  needed  from  which  VG ,  is  calculated 
for  stress,  displacement,  and  buckling. 

Consider  the  finite  element  method  Ku  =  P .   Taking  the 
derivative  of  both  sides  with  respect  to  the  design  variable 
Xq  ,   we  have : 


_L_   [Ku]  =   —   (P)  (Eq.  3) 

3x£    =-      Zxl 

Assuming  the  loads  P  are  not  a  function  of  x„ , 

[JL  (K)]  (u)  +  K  —  (u)  =  0  (Eq.  4) 

3x£  *  =  3x^ 

Finally,  we  arrive  at  (u)  =  -K~   C (K)]  u    (Eq.  5) 

3x£   "  dxz         * 

where  K~   is  the  inverse  of  K.   For  efficiency  K    is  not 

actually  calculated. 


IV.   PROGRAM  FEATURES 

A.   INTRODUCTION 

Each  computer  code  has  its  own  special  features  with  which 
the  user  should  be  familiar  if  the  program  is  to  be  used 
effectively.   The  next  few  sections  will  discuss  an  overview 
of  the  code  and  a  typical  problem  that  might  be  solved  using 
the  program.   This  problem  along  with  other  numerical  examples 
will  be  presented  in  detail  with  results  in  Chapter  VI. 

The  FINITE  ELEMENT  code  was  written  to  be  used  as  a  stand 
alone  analysis  program  or  an  analysis  code  that  could  easily 
be  coupled  to  an  optimizer  (of  the  user's  choice)  through 
simple  modifications  to  the  main  driver  program. 

With  user  supplied  area  and  coordinates  coupled  with  input 
control  parameters,  the  analysis  mode  will  calculate  the  weight 
of  the  truss  structure.   Design  variables  may  be  chosen  as 
member  areas,   joint  coordinates,  or  both.   Additionally, 
gradient  information  will  be  calculated  with  respect  to  area 
variables,  coordinate  variables  or  both.   Coupled  to  an  opti- 
mizer the  code  will  optimize  the  weight  of  the  structure  and 
print  the  final  optimization  information  and  gradient  vectors 
of  specified  constraints. 

The  following  example  of  a  25-bar  space  tower  presented 
in  Tables  (I-IV)  shows  some  of  the  options  the  code  contains. 


19 


Table  I 
25 -BAR  TOWER  SAMPLE  OUTPUT 


INPUT  CONTROL  PARAMETERS 


INPUT    PARAMETERS    FOR    STRUCTURAL    ANALYSIS    ANO    DESIGN    ROUTINE, 
25-3AR    SPACP    TPUSS    (TEST    CASE    STP  ESS/DI  SP./3UCKLI.M3/FREU. ) 


"SAO" 


CONTROL    PARAMETERS 


NUMBER    OF    ELEMENTS*    MS  =  25 

TQTAL    NUMBER    CF    JOINTS*    NJ  =  10 

JOINT    CONSTRAINT    VARIABLE,    NCJ  =  10 

NO.    OF    MA^E RIAL    TYPES,     NMT  =  1 

NO.    OF    LOAD    CCNO.,     NLC  =  2 

NO.    CF    EIGENVALUES,    NEIG  =  2 

NO.    PF    «=ig=njvalUES    CALC,    NEIG1=  2 

NO.    OF    FIXED    MASSES,     NFMASS  =  2 

RUCKLING    CONSTRAINT     ID,    NEUBC  =  1 

NO.    OF    DISPL.    CONST., NDSPLC  =  4 

NO.    OF    FRFQ.    CONST.,     NFREQ  =  1 

LUMPED    MASS    OPTION,     LMASS  =  0 

VA"B.    CALC.    CONTROL    ,     IDVCLC  =  2 


ACCELERATION    OUE    TO    GRAVITY,     GRAV 
EIGENVALUE    CONVERGENCE    TOLORANCE, 


=    0.38640E+03 
EPS6IG   =    0.10000E-Q3 


JOINT    COORDINATES 

JOINT  X 

1  -0.375G0E+02 

2  0.375005*02 

3  -0.3  75  00E  +  02 
<*  0.375G0E+02 

5  0.37500E+02 

6  -0.37500E+02 

7  -0.10000E+03 
3  0.10000E+03 
9  0.100005+06 

10  -0.L0000E+3  3 


0.0 
0.0 

0.37500E  +  02 

0.37500F+02 

•0.37500F+02 

•0.  3  75  00F  +  02 

0.10000E+03 

0.100  00E+03 

•0.10000E+03 

■0.10000E+03 


0.20000E+03 

0.20000Et-03 

0.10000E+03 

0.10000E+03 

0.10000E+03 

G.10000E+03 

0.0 

0.0 

0.0 

0.0 


COORDINATE    DESIGN    VARIABLES 
DESIGN   VARIABLE 


JOINT 

1 
2 
3 
'-* 
5 
6 
7 
8 
9 
10 


X 

0.0 
0.0 

•0.10QOE+01 
0.10C0E+01 
0.1000E+01 

•0.10QUF+-01 

•0.1000E  +  01 
0.1000E+01 
0.1000E+01 

•0.1000E+01 


MULTIPLIER 


Y 

0. 
0. 

0. 
0. 

-0. 
-0. 


10005*01 

1000E+0L 

1000E+01 

1000E+-G1 

0. 1000E+-01 

0.1000E+01 

-0.1000E+01 

-0.1000E.+01 


I 

0.0 

0.0 

0.1000E+-01 

0.10G0E  +  01 

0.1000E+01 

0.1000E+01 

0.0 

0.0 

0.0 

0.0 


MATERIAL    PROPERTIES 
CODE  E  RHO 

1      0.1000E+-C8      Q.1000E  +  00 


SIG-MIN 
-O.+0OOE+05 


S  IG-MAX 
0.4Q00E+05 


K-EULER 
0.3927E+02 
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Table    II 

2 5 -BAR   TOWER    SAMPLE    OUTPUT 

JOINT    DISPLACEMENTS    CALCULATED     BY    SYSTEM    ROUTINE    "SAO" 

LOAD   CONDITION  1 

DISPLACEMENTS  DEGREE    OF    FREEDOM 


JOINT 

1 

2 

3 

1 

-0.219G5E- 

-02 

0.38017E+00 

-0.27099E-01 

2 

0.21910E- 

-02 

-0.33017E+00 

-0.27G99E-01 

3 

0.90789E- 

-01 

-0.15964E-01 

-0.68752E-01 

4 

0.91278E- 

•01 

0.17510E-01 

0.3&L00E-01 

5 

-0.90789E- 

-01 

0.15964E-01 

-0.68752E-01 

6 

-0.91278E- 

-01 

-0.1751 LE-01 

0.36100E-01 

7 

0.0 

0.0 

0.0 

3 

0.0 

0.0 

0.0 

9 

3.0 

0.0 

0.0 

LO 

0.0 

0.0 

CO 

LOAD    CONDITION  2 

DISPLACEMENTS  DEGREE    OF    FRcEOCM 

JOINT  1  2                                     3 

1  0.20126E-01  0.3886DE+00      -C.21G23E-01 

2  0.22911E-01  0.38860E+00      -0.32637E-01 

3  0.99531E-03  0.25O51E-G1       -G.95653E-01 

4  0.64733E-02  0.26707E-01       -0.10297E+00 

5  0.81493F-03  Q.24435E-01          0.b2374E-01 

6  0.66536E-02  Q.25192E-01          G.70194E-01 

7  0.0  0,0  0.0 
3  0.0  0.0  0.0 
9  0.0  0.0                                0.0 

10  0.0  0.0                              0.0 

JOINT  COORDINATES 

JOINT  X  Y                                        L 

1  -Q.37500E+02  0.0  0.2Q00GE+G3 

2  0.37500E+02  0.0  0.200GOE+03 

3  -0.37500E+02  0.375 OOE+02            0.10000E+03 

4  0.375OOE+02  0.375Q0E+02             0. 10000E+03 

5  0.37500E+02  -0.375  OOE  +  02            0. 10000E+03 

6  -0.3750OE+D2  -Q.3  7500E+02            Q.10000E+03 

7  -0.10000E+03  0. 100Q0E+03  0.0 
3  O.lOOOOE+03  O.lOOOOE+03  0.0 
9  0.  10000E  +  03  -0.1QQ00E+O3            0.0 

10  -0.10000E  +  03  -0.130  00E+03            0.0 

JOINT  CHORDINATES 

JOINT  X  Y                                        I 

1  -0.37500E+02  0.0  Q.20000E+03 

2  0.37500E+02  0.0  0.200G0E+C3 

3  -0.37500E+02  0.3  7500E  +  02  0. 1 OOOOE+03 

4  0.37500E+D2  0.37500E+02             0 .  100  J0E+03 

5  0.37500EO2  -0.37500E  +  02             0.10000E+G3 

6  -0.375G0E+-02  -0.37500E+02            0.100G0E+G3 

7  -G.10000E+03  0. 10000E+03  J.O 
3  0.  1000GEO3  0.1  JOOQF-j-03  0.0 
9  0.10000E+03  -0.100005+03             0.0 

10  -0.  1000GE  +  03  -0.  10G00E  +  03             0.0 
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Table    III 
2 5 -BAR   TOWER    SAMPLE    OUTPUT 


AREA/COORDINATE    VARIABLES 


THERE    ARE  8      AREA    DESIGN    VARIABLES 

INITIAL    VALUES 

0.200O0E*01    0.20000E+OL    0.20000E+01    0.20000E+01    0.20G00E+01    0  .20000E  +  0 1 

LOWER    BOUNDS 

O.lOOOOE-04    O.LOOOOE-O*    O.lOOOOE-04    O.lOOOOE-04    O.lOOOOE-04    O.lOOOOE-04 

UPPER    BOUNDS 

O.5000OE+02    0.50000E+02    0.50000E+02    0.5000GE+02    Q.50000E+Q2    0.50000E+02 

THERE  ARE    5   COORDINATE  VARIABLES 

INITIAL    VALUES 

0.37500E+02    0.37500E+02    0. 10000E+03    0.10000E+03    0.L0000E+03 

LOWER    BOUNDS 

O.IOOOOE+OO    O.IOOOOE+OO    O.IOOOOE+OO    O.IOOOOE+OO    O.IOOQOE+OO 

UPPER    BOUNDS 

0.50000E+03    0.50000E+03    0.50OOQE+03   0.50000E+O3    0.50000E+03 


JOINT    DISPLACEMENT    CONSTRAINTS 
DIRECTION    1=X,    2=Y,    3=Z,    0=RESULTANT 


LOAD 

LOWER 

UPPER 

NOOE 

OIR. 

COND. 

BOUND 

BOUND 

I 

1 

1 

-0.35TOE+00 

0.350CE+O0 

2 

1 

I 

-0.3500E+00 

0.350CE+GO 

1 

2 

L 

-0.3500E+00 

0.3500E+30 

2 

2 

I 

-0.3500E  +  00 

0.3500E+00 

FREQUENCY    CONSTRAINTS 
LOWER    BOUNDS     IN   CPS 

FREQUENCY  LOWER 

NUMBER  BOUNO 

L  0.1600E+02 
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Table  IV 
2 5 -BAR  TOWER  SAMPLE  OUTPUT 

FINAL  OPTIMIZATION  INFORMATION 


PINAL    OPTIMIZATION     INFORMATION 


03J 


0.2672566*03 


DECISION    VARIABLES     IX-VECTQ?) 

1)  0.145956*01       0.74285E*00      0.983926*00 

7)  0.85953F*)0       0.69571E  +  00       0.39367E+02 

13)  C.68280E+02 


CONSTRA 
I) 
7) 

13) 

19) 

25) 

31) 

37) 

43) 

'♦9) 

55  ) 

61) 

67) 

73) 

79) 

85  ) 

91  ) 

97) 
103) 
109) 
115) 
121) 
127) 
133) 
139) 
145) 
151  ) 
157) 


INT  VALUE 
-0.238b6E 
-0.29252E 
-0.10090E 
-0.39607F 
-0.91932E 
-0.10259E 
-0.13491E 
-0.75447E 
-0.11021= 
-0.766  VOF 
-0.U145E 
-0.10366F 
-0.13589E 
-O.U013E 
-0.923646 
-0.1  I726c 
-0.806f7E 
-0.11518E 
-0.78602E 

-c.840oos 

-0.8U20E 
-0.U316E 
-C.11028E 
-0.12234F 
-0.72096E 
-0.68067E 
-0.11831E 


S  (G- 
-01  - 
-01    - 

♦  01  - 
•■00  - 
+  )0    - 

♦  31    - 

♦  01    - 

♦  00    • 

♦  01    - 

♦  JO    - 

♦  01    - 

♦  01    - 

♦  01    - 

♦  01    - 

♦  00  - 
♦01   - 

♦  00   - 

♦  01    - 

♦  00  - 
*-0  0    • 

♦  00    - 

♦  01    - 

♦  01    • 

♦  01  - 
♦00   - 

♦  JO    ■ 

♦  01    - 


vEcn 

0.993 
0.292 
0.991 
0.  110 
■0.U8 
■0.9  74 
0.950 
0.124 
0.397 
0.123 
0.385 
0.963 
0.941 
0.898 
0.107 
0.327 
0.  119 
•0.848 
0.121 
■0.  116 
•0. I  18 
0.36U 
0.897 
0.776 
•0.  127 
•0.131 
•0.316 


P  ) 

17F+00 

52F-01 

OOF+JO 

39F*0l 

07E+Q1 

13E*00 

88E*00 

55E*01 

396*00 

31E+01 

46E*00 

43E*00 

11E+00 

176*00 

64E*0l 

40F*00 

32F*0l 

15E+00 

40E+01 

OOE+Ol 

88E+01 

40E+00 

21E+00 

bie+OO 

906+01 

9  36*01 

9GE*00 


-0.  100 
-0.197 
-0.  100 
•0.631 
•0. 8o9 
•0.104 
•0.  107 
-O.o49 
-0.  106 
-0. 056 
-0.107 
-0.1u3 
-0.  105 
-0.  105 
-0. 957 
-0.  130 
-0. 969 
•0.1  10 
-0.  173 
-0.383 
-0.272 
-0. 150 
-0.139 
-0.  168 
-0.l<t4 
-0.212 
-0.156 


68E*0l 
07E*01 
356*01 
31 E+00 
05E*00 
20E*0l 
97E*01 
20E*00 
27E+01 
83E*00 
036*01 
156*01 
076*01 
70E*01 
236*00 
66E*01 
«9F-G1 
9oE  +  01 
45E-04 
256*00 
196*00 
73E*01 
62E+01 
476*01 
73 E+00 
37E-01 
12E+01 


THERE    ARE  7    ACTIVE    CONSTRAINTS 

CONSTRAINT    NUMBERS     ARE 

I  7  8       111       120       126       153 

THERE    APE  0    VIOLATED    CONSTRAINTS 

THERE    ARE  0    ACTIVE    SIDE    CONSTRAINTS 

TERMINATION    CRITERION 
ITFR    EQUALS     ITMAX 

NUMBER    OF    ITERATIONS    =         20 

OBJECTIVE    FUNCTION    WAS    EVALUATED 

CONSTRAINT    FUNCTIONS    WERE    EVALUATED 

GRAOIENT    CF    OBJECTIVE    WAS    CALCULATED 

GRADIENTS    OF   CONSTRAINTS    WERE    CALCULATED 


0.10425E+01 
0.469346*02 


-0.10068E*01 
-0.1Q160E*01 
-0.71-*o5E+00 
-O.12349E*0l 
-0.  123496*01 
-0.71*656+00 
-0.12113E+01 
-0.  o55366+C0 
-O.o55366*00 
-0. 12113E+G1 

_n      i   r\  j.  7  j  c  j. 


00    - 


-u.  /4a/  /t*uu 
-0.  11951E+01 
-0. 994166*00 
-0.90920E+00 
-0.994166+00 
-0.  90919E  +  00 


60 

TIMES 

60 

TIMES 

20 

TIMES 

20 

TIMES 

0.112716*01 

0. 

735676*00 

0.133956*03 

0. 

/3636E+02 

0.993166*03 

-0. 

19707E*0i 

0.984056*00 

-0. 

10063E+01 

0.  128546*01 

-0. 

53bd4E*00 

0.765076*00 

-0. 

1  38136*01 

0. 765076*00 

-0. 

138136*01 

0.123546*01 

-0. 

52b>i*c  +  OQ 

0.788706*00 

-0, 

,11298E*J1 

0. 13<«-4oE*01 

-0. 

7b832E*03 

0.134466*01 

-0. 

,  7bd32E*00 

0.  786706*00 

-0. 

112986*01 

0. 952166  +  00 

-0. 

104126*01 

0.9521*6*00 

-0. 

,104126*01 

0.98304E+00 

-0. 

10095E*0l 

0.983046*00 

-0. 

■10095E+01 

0. 399026+00 

-0. 

,14  7196*01 

0.114236*01 

-0. 

332846*00 

■0.li42dE*01 

-0. 

332o4E*00 

O.d99326*00 

-3. 

lt719E+0l 

0.30493E*00 

-0. 

175206*01 

0.125126*01 

-0. 

,31d50E-01 

0.125126*01 

-0. 

,315516-01 

0.80493E*03 

-0. 

l/52CE*0l 

0.10053E*01 

-0, 

,982096*00 

0.10908E+01 

-0. 

721686+00 

0.100586*01 

-0. 

982096*00 

0.109086*01 

-0. 

721S8E+00 
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V.       USER   GUIDE 

A.  INTRODUCTION 

In  developing  any  computer  code  for  engineering  analysis, 
it  is  necessary  to  additionally  develop  concise,  easily  under- 
stood software.   This  USER  GUIDE  is  writen  to  be  easily 
followed  assuming  minimal  FORTRAN  knowledge.   The  format 
follows  that  of  the  optimization  code,  COPES/CONMIN , 
reference  1. 

This  chapter  is  devoted  to  acquainting  the  user  with  the 
code  and  necessary  input  data.   A  simple  3-bar  truss  analysis 
is  used  as  the  example . 

B.  GENERAL  FORMULATION 

Each  data  card  is  set  up  to  show  the  input  data  necessary 
with  the  3-bar  truss  example  underneath.   Formats  are  of  "15" 
and  "F10.0"  type.   "I"  formats  must  be  right  justified,  and 
"F"  formats  must  have  the  decimal  point.   The  number  of  cards 
read  per  data  block  is  listed  at  the  bottom  of  each  block. 

C.  CONSTRAINTS 

Constraints  are  calculated  and  stored  in  the  G  vector  as 
listed  in  the  following  chart.   The  total  number  of  constraints 
NCON  =  2-NDSPLC  +  NFREQ  +  NE*NLC .   When  any  of  the  constraints 
are  missing  from  the  G  vector,  all  constraints  are  moved  up. 
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For  example,  if  there  is  no  frequency  constraint,  then  a  diS' 
placement  constraint  would  fill  the  first  location  of  the  G 
vector . 


G  = 


-  1 


NFREQ 

(FREQUENCY  CONSTRAINTS) 


2*  NDSPLC  *  NLC 
(DISPLACEMENT  CONSTRAINTS) 


MEMBER  #1 
LC  #1 


NE 


_  i 


MEMBER  #2 
LC  #1 


MEMBER  #NE 
LC      #NLC 


STRESS  AND 
BUCKLING 


STRESS  AND 
BUCKLING 


STRESS  AND 
BUCKLING 
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D.   EXAMPLE 

The  initial  layout  of  the  3-bar  truss  is  shown  in  Figure 
6.1.   Stress  constraints  were  imposed  as  well  as  constraints 
on  Euler  buckling,  displacement,  and  first  fundamental 
frequency . 

With  geometry  specified  as  per  the  figure,  2  independent 
load  conditions  (P1,P2)  were  imposed  with  3  member  sizing 
variables  (Al,  A2 ,  A3)  linked  so  that  A1=X(1),  A2=X(2), 
A3=X(1) . 

1.  Properties/Conditions 

3 
The  material  used  had  a  density  (p)=0.1  lb/in  . 

7 
Young's  Modulus  was  selected  as  10   psi. 

LOADS:  P1=P2=  20,000  lbs. 

STRESS  LIMITS:  -15,000  <  o.  <  20,000  psi     i=l,3 
Symmetry  was  maintained  by  linking  of  variables  so 
A1=A3 .   Initial  design  began  with  A1=A2=A3=  1.0  in.   One  non- 
structural mass  of  500  lbs  was  attached  at  joint  2  .   Displace- 
ment constraints  were  imposed  between  -.02  and  .02  inches  in 
the  x-direction  with  the  first  natural  frequency  limited  to 

a  value  u>   <  20  Hz. 
n 

2 .  Input  Control  Parameters 

The  following  input  control  parameters  are  given  for 
ease  of  the  following  example.   NE=3,  NJ=4,  MCJ=4,  NMT=1,  NLC=2, 
NEIG=2,  NEIG1=2,  NFMASS=1,  NEUBC=1,  NDSPLC=2,  MFREQ=1,  LMASS=0, 
IDVCLC=2.   Table  V  is  a  listing  of  commonly  used  nomenclature. 
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TABLE    V 

COMMON   VARIABLE    NOMENCLATURE 


A-  MEMBER1  S  CROSS-SECTIONAL  AREA 
BL-  LOWER  BOUND  ON  DISPLACEMENTS 
BU-  UPPER  BOUND  ON  DISPLACEMENTS 
DIR-    DIRECTION    1=X,       2=Y,        3=Z. 

0=LIMIT    ON    RESULTANT    DISPL. 
E-    YOUNGS    MODULUS 
ELNO-    ELEMENT    NUMBER 
EPSEIG-    CONVERGENCE    TOLERANCE   OF    EIGENVALUE 


SOLUTIONS     (DEFAULT  =  .0001) 

, FCN-    LOWER    BOUND   ON    Pl«  .   .  ,      _ 

SYSTEM    NATURAL    FREQUENCY 


FC1.FC2,.  .  ..FCN-    LOWER    BOUND   ON    FIRST,  SECOND ,    ETC 


FX-    LOAD    FORCES    BEING   APPLIED    IN    THE    X    DIRECTION 
FY-    LOAD    FORCES    BEING   APPLIED    IN    THE     Y    DIRECTION 
FZ-    LOAD    FORCES    BEING   APPLIED    IN    THE    Z    DIRECTION 
GRAV-    ACCELERATION    DUE    TO    GRAVITY     (DEFAULT=    386.4) 
IX-    CONSTRAINT    IDENTIFIER.     IF    NON-ZERO    THE    X-DOF 

IS    CONSTRAINED 
IY-    CONSTRAINT    IDENTIFIER.     IF   NON-ZERO    THE    Y-DOF 

IS    CONSTRAINED 
IZ-    CONSTRAINT    IDENTIFIER.     IF    NON-ZERO    THE    Z-DOF 

IS    CONSTRAINED 
JN-    JOINT    NUMBER     (GLOBAL) 
KEULER-    EULER    BUCKLING    COEFFICIENT    FOR    BAR 

ELEMENTS 
LC-    LOAD    CONDITION 
LMASS-    LUMPED    MASS    OPTIONS.    IF    LMASS    . NE.     0    LUMPED 

MASS    MATRIX    USED 

IF    LMASS    =0    CONSISTENT    MASS    MATRIX    USED 
NEIG1-NUMBER    OF    EIGENVALUE/EIGENVECTOR    TO    BE 
EVALUATED    DEFAULT    =  MIN(2*NEIG,     NEIG+8) 
NEUBC-    BUCKLING    CONSTRAINT    IDENTIFIERS     .     IF    NEUBC. 

NE.O    EULER    BUCKLING     WILL    BE    IMPOSED    IN 

EAR    ELEMENTS 
NFMASS-NUMBER    OF    FIXED    MASSES    ATTACHED    TO    THE 

STRUCTURE 
NFREQ-NUMBER    OF    FREQUENCY    CONSTRAINTS 
NID-    NUMBER    OF    INDEPENDENT     DEGREES    OF    FREEDOM 
NJ-NUMBER    CF    JOINTS 
NLC-    NUMBER    OF    LOADING   CONDTIONS 

NLJ-    NUMBER    OF    LOADED   JOINTS    FOR    THIS    LOAD    CONDITION 
NMT-    NUMBER    OF    SEPARATE    MATERIAL    TYPES 
RHO-    MATERIAL    DENSITY 
SIGMIN-    MINIMUM    ALLOWABLE    STRESS 
SIGMAX-    MAXIMUM    ALLOWABLE     STRESS 
XA-    INITIAL    VALUE    OF    AREA    DESIGN    VARIABLE 
XAL-    LOWER    BOUND 
XAU-    UPPER    BOUND 

XC-    INITIAL    VALUE    OF    COORDINATE    DESIGN    VARIABLE 
XCL-    LOWER    BOUNDS 
XCU-    UPPER    BOUNDS 
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The  following  user's  manual  is  divided  into  blocks  A 
through  M.   Appearing  directly  below  each  data  field  line 
are  the  parameters  for  the  3-bar  truss  example.   It  is  impor- 
tant to  note  the  user  may  choose  any  units;  however,  all  units 
must  remain  consistent  throughout  the  problem. 
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r 

i 

|       DATA    BLOCK       A 
1 

DESCRIPTION;    Title    Card 

Format   and  Example 

TITLE                                                                               | FORMAT 

1 

1                                                                                                      J20A4 

1 

1 
J 

3-BAR    TRUSS     (EXAMPLE) 

1 

1 

1 

FIELD                                           CONTENTS 

1                                  ANY    80    CHARACTER    TITLE 

MAY    BE    GIVEN    ON    THIS    LINE 

i 
■ 

J 
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DATA  BLOCK   B 


DESCRIPTION:  Control  Parameters 


Format  and  Example 


NJJ  NCjT  NMT|  NLCJNEUBCl  NDSPLCJ  NPRI  |FORMAT 

iii         i  I         L_. 

ii  i 

—i i,       i i 1 


NE 


x. 


1       J    2      j       12 

I I J J_ 


0 


NFREQ    ilDVCLCl  IRECP       1  FORMAT 

1  |       2        |       0  | 515 

I                 I                    I 
j i i 


NEIG 
2 


NEIG1 
2 


JL 


NFMASS 
1 


JL 


T 


I 

I 

-I 


NOTE:    DEFINITIONS    OF     INPUT    CONTROL    PARAMETERS    FOR 
PROGRAM    NEXT    PAGE. 
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1 

HELD 

CONTENT 

1 

NE-number   of    elements 

2 

NJ-number   of    joints 

3 

NCJ-number  of  constrained    joints 

4 

NMT-number   of  seperate    material   types 

5 

NLC-number  of   load  conditions 

6 

NEUBC- buckling  constraint   identifier 

(If   NEUBC.NE.O   -EULER    buckling  con- 

straints will   be    imposed   on    bar    ele.) 

7 

NDSPLC-number   of    displacement   constraints 

8 

NPR1-input   print   control    (if   NPRI.ne.O 

input  information    will   not   be    printed) 

1 

NEIG-number    of  precise    eigenvalues   to 

be  evaluated 

2 

NEIG1-number    of  eigenvalues   to   be    evaluated 

DEFAULT=min.    of     (2*NEIG    ,    NNEIG  +  8) 

3 

NFMASS-number  of    fixed    masses  attached   to 

structure 

4 

NFREQ- number    of  frequency   constraints 

5 

IDVCLC-design   variable    control   parameter 

If     (IDVCLC.EQ.1)     NDV=NDVAR1 

If     (IDVCLC.EQ.2)     NDV=NDVAR1+    NDVAR2 

If     (IDVCLC.EQ.3)     NDV=NDVAR2 

6 

IRECP-recripocal    variable   identifier 

If    (IRECP.GT.O)    the    X-vector    contains 

the   recripocal   of   area   and   VLB   &    VOB 

are  changed   accordingly 

—  . 
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DATA    BLOCK       C 


DESCRIPTION:    Dynamic   Analysis   Information 


Format   and   Example 


|     FORMAT 


LMASS 


1 


GRAV 


EFSEIG 


j_ 


I 


15, 
2F10.0 


386. U    | 
I  I 


0. 


FIELD 


CONTENTS 

LMASS -lumped   mass  option    (if   LMASS. NE.O) 
the   lumped   mass    matrix   is   used. 

GRAV-accleration   due   to   gravity 
(default=386.4    inches/sec  ) 

EPSEIG-conv erge nee    tolerance   on   eigenvalue 
solution,     (defaults. 0001) 
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DATA  BLOCK   D 

DESCRIPTION:  Joint  Coordinates  S 

Design  Variable  Linking 

Format  and  Example 


jn| 


ix|iyTiz 


PCX 


PCY       I    PCZ    ] FORMAT 
|  Il5, 


i.      i i      i 


3F10.0, 

315, 

j  3F10.0 
I 


1 j 

-10. 



0.     I     0.     I    1)    Oj    Ol    -1. 

1                            III! 

0.     j       0.     j 

L          X-                 J 

-10. 

o.       oj   o]   oj     0.   ]     o.    t     0. 

1         1          1         1                   1                    1                   1 

I3 

• 

0. 

0. 

o.        oj   o]    oj      0.    J      0. 

till                                 ! 

•• 

• 

'    a 

10. 

„. I 

0. 

L—         — I 

o.  TT 

1  ■ 

0       0]       1.1       0.     1       0. 

X X J J— I 

FIELD 


1 
2 
3 
U 
5 
6 
7 
8 
9 
10 


I 


CONTENTS 

JN-joint  coordinate   number 

X-x   coordinate 

Y-y   coordinate 

Z-2   coordinate 

IX-design  variable   associated   with   x   coord. 

IY-design  variable   associated   with   y    coord. 

IZ-design  variable   associated   with   z    coord. 

PCX-participati en   coefficient   of    x-coord. 

PCY-partici pati on   coefficient   of    y-coord. 

PCZ-participati on   coefficient   of    z-coord. 


NOTE:    Number   of   cards  read=NJ 
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DATA    BLOCK       E 


DESCRIPTION:    Material   Properties 


Format   and   Example 


E 

RHO 

SIGMIN 

SIGMAX 

KEULER    |     FORMAT 
j     5F10.0 

I 
I 

i 

1.                                 —J 

L    „                                  1 

L J 

L.  . J 1 

I 

i 

[1000000. i       .1 

l L__ 


I     -1500  0.         20000.            4. 
-i i i j 


FIELD 


1 
2 
3 

a 
5 


CONTENTS 

E-Young's   Modulus 
RHO- material   density 
SIGM IN-minimum    allowable   stress 
SIGM AX-maximum    allowable  stress 
KEUL ER-Euler   buckling   coefficient 
for   bar    element 


NOTE:    Number   of  cards   read=NMT 
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DATA    BLOCK       F 


DESCRIPTION:    Bar   Element   Information 


Format   and   Example 


ELNO    | 
■                   i 

NODE  1J  NODE  2    iMArCODJNDSG      T    A 
iii 

| FORMAT 
■ 

1          1 

■                   i 

i                   i                   i 

5I5/P10.0 
■                       i 

1 

L...     —     L_ 

1                   1 

.,  i i. ....    j i 

-A—           _     | 

FIELD 


TT 


771 


1 

j.  _ 

L 

JL 

1 

JL 

_         J 

[ 

2  i 

• 

3 1  2 

• 

1 

i 

i2 

i 
i 

i 

L 

3 

! 

4       i       2 

~L..          — 

' 

1 

I 

1 

1 

J 

1.    I 

J 

CONTENTS 

ELNO-element   number 

NODE1-global  number   associated   with 

element  node 
NODE2-global  number   associated   with 

element  node 
MATCOD -material   type   of  this    element 
NDSG-design   variable  number    associated 

with   this   element 
A-member  cross-sectional   area 


NOTE:    Number   of  cards  read=NE 


35 


— i 

DATA    BLOCK      G 

DESCRIPTION:    Joint 

Constraint  Information 

1 

Format 

and 

Example 

1  JN  i 

IX     l 

IY 



IZ 

— i 

FORMAT 

1 

i       I 

1 
1 

! 

415 

1 
i 

j 

.  — . 

._ 

.  _  

fH' 

"■'"T" 

I 

1 

□ 

t_     __L 

"T" 
I 

I 

_   .1... 

0 

1  I 

i 

i 

i 

1 

1 

*- 

1 

□ 

■"■                                      i 

1 

1 

L_    _ 

—  ■■«■■       -         

£!£££ 

CO  I 

ITENTS 

1 

JN-; 

joint 

numb  er 

2 

IX 

x,y,z 

constraint   identifier. 

(if 

non- 

3 

IY 

zero 

the    corresponding   degree 

of 

freedom 

4 

IZ 

is   constrained) 

NOT! 

? . 

Number  c 

)f  car 

ds 

read= 

=  NCJ 
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DATA    BLOCK       H 


OMIT   this   card    if   NLC=0    was  rsad    in   block    B 


DESCRIPTION:    Joint    Loading  Information 


Format   and   Example 


r  nlj  t 


1    FORMAT 

I 

I    15 


1    1    I 
I         I 

I 
i 

• 

!   '  I 

I 

! 

j 

JN    I       FX 

I 
I 


) FORMAT 


FY 


FZ 


1 15, 3F 10.0  | 


I"       2  14140.     i     -14140.  |             0.    ~T~ 

II  I                        t                         II 

i i  i j i i 

\        2       |  -14  140.  J     -14140.  ]             0~!       ]                   | 


FIELD 


CONTENT 

NLJ-number  of   loaded    joints 
for    this   load   condition 
JN- joint   number 
FX 

FY-    Forces   in    the    X,Y,Z    directions 
FZ 
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DATA    BLOCK       I 


OMIT   this   block  if  NFMASS  =  0    was   read   in   block    B 
DESCRIPTION;    Lumped   Mass    Information 


Format   and   Example 


JN         T    MASS    ~T 


I 


i  500.  r 


FORMAT 
I5,F10.0 


FIELD 


CONTENTS 


J N- joint    number 

MASS -concentrated   mass   at    joint    (JN) 
mass   is  in   force   units 
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DATA    BLOCK      J 

OMIT    this    block   if   NDVAR1=0 

DESCRIPTION:    Design  Variable    Information 

(AREA    Variables) 

Format   and   Example 


r    XA(D 


T    FORMAT 

|    8F10.0 
I  I 


1. 


1. 


.01  .01 


10.  10. 


u 


j- 


XAL(I) 


1    XAU(I) 


FIELD 


CONTENTS 

X A- initial  value   of    area  design  variables 

XAL-lower   bounds   on    area  design  variables 

XAU-upper    bounds   on    area  design  variables 


NOTE:    read   one   value  of    XAr    XAL,    XAU   for    each 

independent    area    variable  defined    in   Block    D 

Number   of   cards    read    =as   required 
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DATA    BLOCK       K 

OMIT    this    block"  if   NDVAR2  =  0 

DESCRIPTION:    Design   Variable    Information 

(COORDINATE   Variables) 


Format   and   Example 


r         •    - — ' 

XC(1) 

I                                 I8F10.0 
ii                                 ii 

i         10. 

i 

i               !                     I 
i               i                     i 

i            1.0 

I 

1  XCL 

1                        1                                  1 
1                        ■                                  f 

I             20. 

1                                   XCO 

1                     I 

i    ..            4 1 

FIELD  CONTENTS 

1  XC -initial    value   of   coord,    design   variables 

2  XCL- lower    bounds   on   coord,    design   variables 

3  XCO-upper    bounds  on   coord,    design   variables 

NOTE:    read    one    value    of    XC,    XCL,    XCU    for 

each    independent    coordinate   variable 
defined   in  Block    D 

Number  of   cards   read    =as   required. 
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DATA    BLOCK      L 


OMIT   this   block  if  NDSPLC=0   in    Block   A 
DESCRIPTION;    Joint    Displacement   Constraint   Information 


Format   and   Example 


BL  B0"T 


I  FORMAT 
|3I5,2F10.6 


1 


2    ]     1     I   1    ~T  -.02T  .02  i 

[        I      J 


:        2      T    2,1  -.02T    .02~~ 

I I        i L L__L 


FIELD 


1 
2 

3 

4 


CONTENTS 

JN-joint    number 
DIR-direction    1=X    2=Y   3=z 

0=limit  on   resultant  displ. 
LC-load  condition 
BL -lower  bound   on    displacement 

(If   DIR=0   read    0    here) 
3U-upper   bound   on    displacement 


Number   of   cards  read=   NDSPLC 
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DATA    BLOCK       H 


OMIT   this   block  if   NFREQ=0    was   read    in   Block   A 
DESCRIPTION:    Frequency  Constraint   Information 


Format   and  Example 


r    FC1     |    FC2 

LJ_ 


|  FOR 


FC3 


poT] 


1_L 


FCN 


I 


MAT 


8F10.0  I 


FIELD 


N 


CONTENTS 

FC1-    lower    bound   on   first   natural   frequency 
constraint  in   CPS.     (cycles    per    second) 

FCN-    lower    bound  on   NFREQ-th   natural 
frequency   constraint   in  CPS. 
(cycles  per   second) 


NOTE:    OMIT   this   block   if   NFREQ=0   was   read   in    data 
block   B. 

Number  of  cards  read=   as   required 
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VI.   NUMERICAL  EXAMPLES 

A.  INTRODUCTION 

Design  of  planar  trusses  and  3-dimensional  space  towers 
are  presented  here  and  the  corresponding  numerical  results 
are  summarized  to  demonstrate  the  use  of  the  code. 

In  the  examples  given  here,  the  design  variables  are 
member  cross-sectional  areas  and  joint  coordinates.   It 
should  be  emphasized  that  in  practical  design,  the  reciprocal 
of  the  member  areas  is  usually  a  better  choice  for  the  design 
variables.   However,  the  purpose  here  is  to  knowingly   create 
difficult  optimization  problems,  thus  the  choice  of  variables 

The   examples   begin  with   the    3-bar  truss . 

B.  CASE  1:   3-BAR  PLANAR  TRUSS 

The  simple  3-bar  planar  truss,  as  shown  in  Figure  5.1  has 
been  previously  used  for  the  user  guide  example.   This  struc- 
ture was  designed  for  optimum  geometry  subject  to  a  set  of 
two  load  conditions,  buckling  constraints,  displacement  con- 
straints, and  a  lower  bound  on  the  system's  first  natural 
frequency.   The  allowable  stresses  specified  are 

-15000.  <  o.    <    20000.  psi  1=1.,  3 

The  member  areas  were  linked  in  the  following  groups:   A1=A3; 
A2  . 


43 


1 .  Case  la 

Only  the  stress  constraint  was  imposed  for  this  case . 
Final  design  information  is  given  in  Table  VI.  The  number  of 
analysis  for  this  design  was  33. 

2 .  Case  lb 

This  case  includes  the  previous  constraint  plus  con- 
straints on  buckling  (KEULER=4 . 0 ) ,  displacement  (-.02  -.02  in.), 
and  first  natural  frequency  (oo  =20.0).   Results  are  given  in 
Table  VII.   The  number  of  analysis  required  for  this  design 
was  16  with  1.73  seconds  of  CPU  time. 


14140  lbs, 


14140  lbs. 


Figure  5.1   3-BAR  TRUSS 
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C.   CASE  2:   18-BAR  PLANAR  TRUSS 

A  cantiliver  truss,  as  shown  in  Figure  6.2,  has  been  pre- 
viously used  as  a  standard  test  case  for  structural  design 
[Ref.  23.   The  structure  was  analyzed  for  a  single  set  of  load 
conditions  with  allowable  stresses  being 

-20000.  <  a.    <    20000.  psi  i=l,18 

n 

Young's  modulus  was  taken  as  10   psi  with  a  material  density 
P  =  0  .1  lb./cu  in. 


©       ' 


© 


250" 


Figure  6.2   18-BAR  TRUSS 
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The  independent  coordinate  variables  were  taken  as  X3,  Y3, 
X5,  Y5,  X7,  Y7,  X9 ,  Y9 .   The  member  areas  were  linked  as 
follows:   A1=A4=A8=A12=A16;   A2 =  A6 =A10 =A14=A18 ;   A3=A7=A11=A15 ; 
A5=A9 =A13=A17 .   There  are  a  total  of  four  independent  area 
variables  and  eight  coordinate  variables . 

1.   Case  2a 

This  design  analysis  includes  stress,  buckling  (KEULER= 
4.0),  displacement  (-10.0  to  10.0  in.),  and  a  first  fundamental 
frequency  of  3  Hz.   Additionally,  a  nonstructural  mass  of  W=5000 
lbs.  was  attached  to  node  1.   The  number  of  analyses  for  this 
design  was  62.   Results  are  presented  in  Table  VIII.   26.72 
seconds  of  CPU  time  were  used  for  this  calculation. 

D.   CASE  3:   2 5 -BAR  SPACE  TOWER 

© 


Figure  6.3   25-BAR  SPACE  TOWER 
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The  25-bar  space  tower  shown  in  Figure  6.3  was  designed 
for  two  independent  load  conditions  given  in  Table  X.   The 
allowable  stresses  were  specified  as  follows 

-40000.  <  a.  <  40000.  psi  i=l,25 

7 
Young's  modulus  was  selected  as  10   psi  with  a  material  density 

p=.l  lb/cu.  in.   Members  are  assumed  tubular  with  a  nominal 
diameter  to  thickness  ratio  of  D/t=100  giving  KEULER= 39 . 2 74 . 
Symmetry  was  imposed  in  both  the  x-z  and  y-z  planes.   Non- 
structural masses  of  W=500  lbs.  were  attached  at  nodes  1  and 
2.   Coordinate  variables  were  X4  ,  Y4,  Z4,  X8,  and  Y3  with  the 
remaining  coordinates  linked  to  maintain  symmetry.   Area  vari- 
ables were  linked  in  the  following  manner:   Al ;  A2=A3=A4=A5; 
A6=A7=A8=A9;  A10=A11;  A12=A13;  A14 =A15 =A16 =A17 ;  A18=A19 =A20 =A21 ; 
A22=A23=A24=A25 . 
1 .   Case  3a 

Stress,  displacement  (-0.35  to  0.35  in.),  Euler  buck- 
ling and  first  natural  frequency  limited  to  a  value  w   >  16 
Hz.  were  imposed  for  this  case.   Final  design  information  is 
given  in  Table  XI.   The  number  of  analyses  for  this  design  was 
61  with  2  7.99  seconds  of  CPU  time  used. 

E.   CASE  4:   2  34-BAR  SPACE  TOWER 

The  inital  layout  of  the  tower  is  shown  in  Figure  6.4 
stress  limits  were  as  follows 

-15000.  <  a.  <  20000.  psi  i=l,234 

7 
Young's  modulus  was  chosen  as  steel  3.xl0   psi  with  a  material 

density  of  aluminum  p  =.1  lbs/cu.  in. 
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1.   Case  4a 

Stress  constraints  as  well  as  constraints  on  Euler 
buckling  were  imposed.   The  resulting  problem  had  56  area 
design  variables  and  42  coordinate  variables.   This  problem 
required  91  structural  analyses  using  72  minutes  of  CPU  time. 
The  objective  function  was  evaluated  89  times  and  gradients 
were  calculated  30  times.   Results  are  presented  in  Tables  XII 
through  XIV.   Although  the  weight  of  the  structure  was  only 
reduced  by  2%,  it  should  be  noted  that  constraints  were  ini- 
tially violated  and  the  optimizer  overcame  the  constraint 
violations .   It  is  believed  the  optimum  that  was  reached  is 
not  the  true  optimum  because  of  the  extreme  non-linearity  of 
the  problem. 
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FRONT/BACK 


LEFT/RIGHT 


Figure  6.4   234-BAR  SPACE  TOWER 
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BODY 

LOAD 
COND. 

1 

2 


Table  VI 
3-BAR  TRUSS.   DESIGN  INFORMATION  (STRESS) 


JOINT 

2 
2 


X 

0  .1414E  +  05 
•0  .1414E  +  05 


LOADS 

Y  Z 

-0.1414E+05  0.0 

-0.1414E+05  0.0 


FIXED  MASSES 

NO .  JOINT 

1  2 


MASS 
500  lbs 


MEMBER 
A1  =  A3 
A2 


AREA  (sq.  in.) 
INITIAL 
1.0 
1.0 


FINAL 

.7662E+00 

.4833E+00 


JOINT 
2 


COORDINATES  (in.) 
INITIAL 
10 


FINAL 
.9856E+01 


FINAL  WEIGHT  FOR  AREA  AND  COORDINATES 
WEIGHT  =  2.6  349  (lbs.) 
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BODY 


Table  VII 

3-BAR  TRUSS.   DESIGN  INFORMATION 
(STRESS,  DISPLACEMENT,  BUCKLING,  FREQUENCY) 

JOINT  LOADING  INFORMATION 
SAME  AS  TABLE  VI 

JOINT  DISPLACEMENT  INFORMATION 


JOINT 

DIR. 

LOAD 

LOWER 

UPPER 

NUMBER 

COND. 

BOUND 

BOUND 

2 

1 

1 

-.02 

.02 

2 

2 

1 
AREA  (sq.  in 

-.02 
.) 

.02 

MEMBER 

INITIAL 

FINAL 

A1  =  A3 

1.0 

0 

.99995E+00 

A2 

1.0 
COORDINATES  ( 

in.) 

0 

.1000  E-01 

JOINT 

INITIAL 

FINAL 

2 

10  .0 

0 

.99986E+01 

FINAL  WEIGHT  =  2.8  3  81  (lbs.) 
CRITICAL  CONSTRAINT  -  FREQUENCY 
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BODY 


Table  VIII 
18-BAR  TRUSS.   LOAD  CONDITIONS 


LOAD  CONDITIONS  (lbs.) 


LOAD  COND. 

JOINT 

X 

LOADS 
Y 

Z 

1 

1 

0  .0 

-0  .2000E  +  05 

0  .0 

1 

2 

0  .0 

-0  .2000E  +  05 

0  .0 

1 

4 

0  .0 

-0 .2000E+05 

0  .0 

1 

6 

0  .0 

-0  .2000E+05 

0  .0 

1 

8 

0  .0 

-0  .2000E  +  05 

0  .0 

FIXED  MASS  INFORMATION 
NO.  JOINT  MASS 

1  1  5000  (lbs.) 
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Table  IX 

18-BAR  TRUSS.   DESIGN  INFORMATION 
(STRESS,  DISPLACEMENT,  EULER  BUCKLING,  FREQUENCY) 


AREA  (sq .  in  .  ) 


MEMBER 

INITIAL 

A1=A4=A8=A12=A16 

10  .0 

A2=A6=A10=A14=A18 

10  .0 

A3=A7=A11=A15 

10  .0 

A4=A9=A13=A17 

10  .0 

FINAL 
38  .20 
38.27 
15  .98 
20  .62 


COORDINATES  (in.) 


JOINT 

3 

5 
7 
9 


INITIAL 
X 

Y 

FINAL 
X 

Y 

1000  . 

0. 

859  . 

0 

750  . 

0  . 

628. 

0 

500  . 

0  . 

384. 

0 

250  . 

0. 

148. 

0 

FINAL  WEIGHT  =  13,940  (lbs.) 
CRITICAL  CONSTRAINT  -  FREQUENCY 
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BODY 


NO. 
2 


Table  X 
25-BAR  TRUSS.   LOAD  CONDITIONS 


LOAD  CONDITIONS  (lbs.) 


LOAD  COMD. 

JOINT 

X 

LOADS 
Y 

z 

1 

1 

0  .0 

2000  . 

-5000 

2 

0  .0 

-2000  . 

-5000 

2 

1 

1000  . 

10000  . 

-5000 

2 

0. 

10000  . 

-5000 

3 

500  . 

0. 

0  . 

6 

500  . 

0. 

0. 

FIXED  MASS  INFORMATION 


JOINT 

MASS 

1 

500 .(lbs  .) 

2 

500  .(lbs  .) 
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Table  XI 

25-BAR  TRUSS.   DESIGN  INFORMATION 
(STRESS,  DISPLACEMENT,  EULER  BUCKLING,  FREQUENCY) 


AREA  I 

'  sq .  in . ) 

MEMBER 

INITIAL 

Al 

2 

.0 

A2=A3=A4=A5 

2 

.0 

A6=A7=A8=A9 

2 

.0 

A10=A11 

2 

.0 

A12=A13 

2 

.0 

A14=A15=A16: 

=A17 

2 

.0 

A18=A19=A20= 

=A21 

2 

.0 

A22=A23=A24: 

--A2  5 

2 

.0 

FINAL 
1.46 
0  .74 
0  .98 
1.04 
1.12 
0  .74 
0  .86 
0  .70 


JOINT 


COORDINATES  (in.) 


X 

INITIAL 
Y 

z 

X 

FINAL 
Y 

Z 

37.5 

37.5 

100 

39  .4 

46  .9 

138.9 

100  .0 

100  .0 

0  .0 

73.6 

68  .3 

0  .0 

FINAL  WEIGHT  =  26  7  (lbs.) 
CRITICAL  CONSTRAINT  -  FREQUENCY 
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Table  XII 
2  34-3AR  SPACE  TOWER.   LOADING  INFORMATION 
BODY 

LOAD  CONDITIONS  (lbs.) 


LOAD  COND. 

JOINT 

X 

LOADS 
Y 

z 

1 

49 

6000  . 

-20000  . 

0  .0 

1 

50 

6000  . 

-20000  . 

0  .0 

2 

49 

6000  . 

-20000  . 

0  .0 

2 

50 

-6000  . 

-20000  . 

0  .0 

3 

49 

6000  . 

-20000  . 

0  .0 

3 

50 

3000  . 

-10000  . 

5000 

4 

49 

3000  . 

-10000  . 

-5000 

4 

50 

3000  . 

-10000 . 

5000 

5 

49 

-3000  . 

10000  . 

5000 

5 

50 

-3000  . 

10000  . 

-5000 

NO 

1 

2 


FIXED  MASS  INFORMATION 


JOINT 

MASS 

49 

200 .(lbs  .  ) 

50 

200. (lbs.) 
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Table  XIII 
2  34-BAR  SPACE  TOWER.   DESIGN  INFORMATION 


AREA  (sq.  in.) 


MEMBER 

INITIAL 

FINAL 

A1,A2 ,A3,A4 

25 

.0 

25  .38 

A5 ,A6 ,A7,A8 

25 

.0 

25  .35 

A9,A10,A11,A12 

25 

.0 

25  .28 

Ak3,A14,A15,A16 

25 

.0 

25  .21 

A17 

,A18,A19,A20 

25 

.0 

25  .14- 

A21. 

A22,A23,A24 

25 

.0 

25  .07 

A2  5. 

,A26 ,A27,A28 

25 

.0 

25  .00 

A29 

,A30,A31,A32 

25 

.0 

24.95 

A33. 

A34,A35 ,A36 

25 

.0 

24  .  89 

A37 

,A38,A39 ,A40 

25 

.0 

24.85 

A41 

,A42 ,A43,A44 

25 

.0 

24  .98 

A45. 

,A46,A47,AU8 

25 

.0 

24  .87 

A49 

,A50 ,A51,A52 

25 

.0 

24.88 

A5  3. 

(A54 ,A55,A56 

25 

.0 

24.89 

A57 

,A58,A59,A60 

25 

.0 

24.89 

A61. 

,A62 ,A63,A54 

25 

.0 

24  .90 

A6  5. 

,A66 ,A67,A68 

25 

.0 

24.91 

A69. 

,A70 ,A71,A72 

25 

.0 

24  .91 

A73 

,A74,A75,A76 

25 

.0 

24  .91 

A77 

,A78,A79 ,A80 

25 

.0 

24  .91 

A81 

,A83 

25 

.0 

24.96 

A82. 

,A8U 

25 

.0 

24.96 

A85 

,A87 

25 

.0 

24.96 

A86 

,A88 

25 

.0 

24.95 

A89. 

,A9  0,A91,A9  2 ,A9  3,A9  4,A9  5 ,A9  6 

25 

.0 

24.68 

A9  7 

,A9  8,A9  9,A10  0 ,A101 ,A10  2 ,A10 3 ,A10  4 

25 

.0 

24  .69 

A150 ,A106 ,A107,A108,A109 ,A110 ,A111 ,A112 

25 

.0 

24.71 

All  3 ,A114 ,A115 ,A116 ,A117 ,A118 ,A119 ,A120 

25 

.0 

24.73 

A121,A12  2,A12  3,A12  4,A12  5 ,A12  6 ,A12  7,A12  8 

25 

.0 

24  .74 

A12  9 ,A130 ,A131,A13  2 ,A13 3 ,A134 , A13  5 ,A136 

25 

.0 

24  .75 

A13  7,A138,A139 ,A140 ,A141 ,A14  2 , A1U 3 ,A144 

25 

.0 

24  .74 

A14  5,A14  6,A14  7,A14  8,A149 ,A15  0 , A151 ,A152 

25 

.0 

24  .72 

A15  3,A15  4,A15  5  ,A15  6  ,A15  7  ,A15  8  ,A15  9  ,A16  0 

25 

.0 

24.69 

A161,A16  2,A16  3,A16  4,A15  5 ,A16  6 ,A16  7,A16  8 

25 

.0 

24.65 

A16  9 ,A170,A171,A172,A17  3,A174,A17  5,A176 

25 

.0 

24.32 

A177,A173,A179 ,A180 

25 

.0 

24  .96 

A181,A182,A183,A184 

25 

.0 

24  .96 

A185,A186,A187,A1S8 

25 

.0 

24  .87 

A18< 

3 ,A190 ,A191,A192 

25 

.0 

24.86 

57 


Table  XIII  (cont'd) 


AREA  (sq .  in . ) 


MEMBER 

A193,A194,A195  ,A196 

A197,A198,A199,A200 

A201,A202 

A203,A204 

A205,A206 ,A207,A208 

A209 ,A210,A211,A212 

A213,A214 

A215,A216 

A217,A218 

A219,A220 

A221,A222 

A223,A224 

A225 ,A226 

A227,A228 

A229 ,A230 

A231,A232 

A233,A234 


INITIAL 

FINAL 

25  .0 

24.95 

25  .0 

24  .97 

25.0 

24  .91 

25  .0 

24  .92 

25  .0 

24.86 

25  .0 

24  .87 

25  .0 

24.89 

25.0 

24.90 

25  .0 

24.91 

25  .0 

24.92 

25.0 

24.92 

25  .0 

24  .93 

25  .0 

24.93 

25  .0 

24  .94 

25  .0 

24  .94 

25  .0 

24.94 

25  .0 

24.93 

COORDINATES  (in.) 


JOINT 

1 

5 

9 
13 
17 
21 
25 
29 
33 


X 

120 
111 
102 
93 
84 
75 
66 
57 
48 


INITIAL 
Y 

0* 

120 
240 
360 
480 
600 
720 
840 
960 


Z 

120 

111 

102 

93 

84 

75 

66 

57 

48 


X 

FINAL 
Y 

Z 

120  .5 

0 

121.4 

110  .7 

120  .4 

112  .1 

102  .2 

241.2 

103.2 

93.1 

362  .6 

93  .7 

84  .0 

483  .8 

84  .4 

74  .9 

604.5 

75  .2 

65  .9 

724.2 

65  .9 

56  .8 

842  .6 

56  .8 

47.8 

960  .4 

47.8 

Initial  weight  =   84,524   lbs. 


Final  weight  =   83,142   lbs, 


*Final   weight    (AREA  VARBIABLES    ONLY)    =    28,900    lbs. 

*  This  clearly  indicates  optimum  not  reached  with  previous   case. 
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Table  XIII  (cont'd) 

COORDINATES  (in.) 


JOINT  INITIAL  FINAL 

X        Y  Z 

37  39  1080  39  38.8  1077.4  38.8 

41  30  0*  30  29.9  -  29.9 


X 

INITIAL 
Y 

Z 

39 

1080 

39 

30 

0* 

30 

30 

1248 

30 

0* 

1248 

90 

15 

1224 

60 

15 

1248 

60 

45  30  1248  30  29.3  1212.0  29.9 

49  0*  1248  90  -  1203.5  90.0 

51  15  1224  60  15.0  1194.1  58.9 

53  15  1248  60  15.0  1214.7  59.2 


-NOT  A  DESIGN  VARIABLE 
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VII.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  finite  element  code  was  presented  coupled  to  an 
optimizer  for  truss  analysis  and  design.   Trusses  were 
designed  for  minimum  weight  with  multiple  load  conditions 
considered . 

The  displacement  method  for  static  analysis  and  the  sub- 
space  iteration  method  for  eigenvalues  were  applied. 

Several  examples  were  considered.   In  every  case  the  code 
worked  as  an  analysis  tool,  and  significant  weight  reductions 
were  obtained  with  the  coupled  optimizer  CONMIN.   Run  times 
and  results  compared  to  the  same  test  cases  with  other  codes 
indiate  that  the  code  is  competitive  as  a  design  tool. 

B.  RECOMMENDATIONS 

The  following  recommendations  may  be  of  value  for  follow 
on  work: 

1.  The  code  should  be  modified  so  the  user  can  access 
stresses  directly. 

2.  The  code  should  be  extended  to  include  other  elements 
such  as  frames  and  plates . 

3.  An  out  of  core  equation  solver  should  be  added. 

4 .  The  method  of  gradient  calculation  should  be  dependent 
on  specific  gradients  required  reference  3  and  reference 
4  . 
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5 .  Gradients  of  frequency  constraints  would  benefit  from 
a  more  efficient  algorithm  reference  5 . 

6 .  The  need  for  a  large  scale  public  structural  optimiza- 
tion code  still  exists  . 
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APPENDIX  A 


DATA  FILES 


A.   INTRODUCTION 

This  appendix  contains  the  data  files  used  to  create  the 
test  cases  in  Chapter  VI.   Additionally  the  data  file  for  the 
user's  guide  in  complete  form  is  presented. 
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(STRESS) 


Table    XIV 
DATA    FILE    3-BAR   TRUSS 


3-    BAR    TRUSS    (TEST    CASE    A-STRESS    ONLY) 


3 

4            4 

1 

2 

0           0 

0 

0 

0            1 

0 

2 

0 

386.4 

0.0 

1 

-10. 

0. 

0. 

1 

0 

0 

-1. 

2 

0. 

-10. 

0. 

0 

0 

0 

0 

3 

0. 

0. 

0. 

0 

0 

0 

0. 

4 

10. 

0. 

0. 

1 

0 

0 

1. 

lOOOOOOO.         .1 

15000         :0OG0. 

4. 

1 

1           2 

1 

1. 

2 

3           2 

2 

1. 

3 

4           2 

1 

1. 

1 

1            1 

- 

2 

0           0 

3 

1            I 

4 

1           1 

1 

2 

14140.      - 

•14140. 

0. 

1 

2 

-14140.       - 

141-40. 

0. 

1 

2 

500. 

1. 

1. 

.01 

.01 

10. 

10. 

10. 

i 

1.0 

20. 
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Table  XV 
DATA  FILE  3-3AR  TRUSS 

(STRESS,  DISPLACEMENT,  BUCKLING,  FREQUENCY) 


3-    BAR 

TRUSS    (TE 

ST  CASE    3- 

S" 

•RE 

SStOISP 

L. ,3u; 

KL.  »P 

RfcQ. 

J 

3 

4           4 

1           2 

1 

2 

0 

2 

2             1 

1           2 

0 

336.4 

0.0 

1 

-10. 

0. 

0. 

1 

0 

0 

-1. 

2 

0. 

-10. 

0. 

0 

0 

0 

0 

3 

0. 

0. 

0. 

0 

0 

0 

0. 

4 

10. 

0. 

0. 

1 

0 

0 

1. 

10000000.         .1 

-15000, 

1 

20000. 

4. 

1 

1           2 

1           1 

1 

. 

2 

3           2 

1           2 

1 

. 

3 

4           2 

1           1 

1 

. 

1 

1            1 

" 

2 

0           0 

3 

1            I 

4 
L 
2 
1 
2  • 

1            1 

14140.      - 

14140. 

0 

• 

-14140.      - 

14140. 

0 

* 

2 

500. 

1. 

1. 

.01 

.01 

10. 

10. 

. 

10. 

' 

1.0 

20. 

2 

1           1 

-.02 

.02 

2 

2           1 

-.02 

.02 

20. 
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Table    XVI 
DATA    FILE    18 -BAR   TRUSS 
(STRESS,    DISPLACEMENT,    BUCKLING,    FREQUENCY) 


18-RA 
13 

2 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
1000000 

1 

2 
3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 

5 

1 

2 

4 

6 

8 

1 
30. 
1. 

1000. 
1000. 
700. 
1200. 

1 
4.0 


R    TRUSS    CASE    (TEST    CASE    STRESS/DISP. /BUC<L ING/FREQ . J 


11  11 
4  1 

386.4 

250. 

000. 

000. 

50. 

50. 

00. 

00. 

50. 

50. 


0.    .1 

1 
1 
2 
2 
3 
3 
4 
4 
5 
5 
6 
6 
7 
7 
8 
8 
9 
9 
0 


250. 
250. 

250. 

0. 

250. 

0. 

250. 

0. 

250. 


1 
1 
0. 


1 


2 
3 

3 
4 
4 
5 
5 
6 
6 
7 
7 
8 
8 
9 
9 
10 
10 

11 

0 


•20000. 
1  10. 
2 
3 
1 
4 
2 
3 
1 
4 
2 
3 


10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 
1  10. 
4    10. 

10. 

10. 

10. 

10. 

10. 


5000. 
30. 
1. 

1000. 
0. 

-50. 

200. 

2  1 


-20000. 
-20000. 
-20000. 
-20000. 
-20000. 


30. 
1. 

1000. 
750. 
500. 
1000. 
-10. 


1 

2 

0    1. 

3 

4 

0    1. 

5 

6 

1. 

7 

8 

1. 

20000. 


4. 


30. 
1. 

1000. 
0. 

-50. 
150. 
10. 


50D. 
200. 
750. 


0. 

-50. 

150. 


25  0. 

50. 
500. 


0. 

-100. 

150. 
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Table  XVII 
DATA  FILE  2 5 -BAR  SPACE  TOWER 
(STRESS,  DISPLACEMENT,  BUCKLING,  FREQUENCY) 


25-8AR 

SPACE 

TRUSS 

(TEST    CASE    STRESS/OISP 

./8UCKLING/FREQ.J 

25 

10 

10 

1 

2 

14            0 

2 

2 

2 

1 

2 

0 

386. 

W 

0. 

1 

-37.5 

0.0 

200.                 0 

0            0 

0 

0 

0. 

2 

37.5 

0.0 

200.                 0 

0            0 

3 

-37. 

37.5 

100.                  1 

2           3 

-1. 

1. 

1. 

4 

37.5 

37.5 

100.                  1 

2           3 

1. 

1. 

1. 

5 

37.5 

-37.5 

100.                 1 

2           3 

1. 

~  L  • 

1. 

6 

-37.5 

-37.5 

100.                 1 

2           3 

"*"  L  • 

~L  ♦ 

1. 

7 

-100 

• 

100. 

0.                      4 

5           0 

™"  X  • 

i. 

0. 

8 

100. 

100. 

0.                      4 

5           0 

i. 

i. 

9 

100 

-100. 

0.                      4 

-      5           0 

i. 

—  1 . 

10 

-100 

-100. 

0.                      4 

5           0 

— i. 

^*  • 

lOOOOOOO. 

".1 

-4000    . 

40000. 

39.274 

1 

1 

2 

1 

2. 

2 

1 

4 

2 

2. 

3 

2 

3 

2 

2. 

4 

1 

5 

2 

2. 

5 

2 

6 

2 

2. 

6 

2 

4 

3 

2. 

7 

2 

5 

3 

2. 

8 

1 

3 

3 

2. 

9 

1 

6 

3 

2. 

10 

3 

6 

4 

2. 

1 

4 

5 

4 

2. 

/ 

4 

3 

5 

2. 

13 

5 

6 

5 

2. 

14 

3 

10 

6 

2. 

15 

6 

7 

6 

2. 

16 

9 

4 

6 

2. 

17 

8 

5 

6 

2. 

13 

7 

4 

7 

2. 

19 

8 

3 

7 

2. 

20 

10 

5 

7 

2. 

21 

9 

6 

7 

2. 

22 

10 

6 

8 

2. 

23 

7 

3 

8 

2. 

24 

8 

4 

8 

2. 

25 

1 
2 

9 

5 

8 

2. 

3 

4 

5 

6 

7 

1 

1 

1 

8 

I 

1 

I 

9 

1 

1 

1 

10 

1 

1 

1 

2 

1 

0. 

20  000. 

5000. 

2 

0. 

-20000. 

5000. 

4 

1 

1000 

10C00. 

5000. 

2 

0. 

10000. 

5000. 

3 

500. 

0. 

0. 

6 

500. 

0. 

0. 

1 

500. 

2 

500. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

0.00001         0. 

00001 

0. 

30001 

0.00001 

0.  00001 

0.00001 

0.00001 

0.00001 

50. 

50. 

50. 

50. 

50. 

50. 

50. 

50 

37. 

5 

37.5 

100. 

100. 

100. 

.1 

.1 

.1 

.1 

.1 

500. 

500. 

500. 

500. 

500. 

1 

1 

1 

-.35 

.35 

2 

1 

1 

-.35 

.35 

1 

2 

1 

-.35 

.35 

, 

2 

2 

1 

-.35 

.35 

_ 

16.0 

.. 
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Table  XVIII 
DATA  FILE  2 34 -BAR  SPACE  TOWER 


234-BAR 

SPACE 

TOWER  DESIGN 

(TEST  CASE 

STRESS»BUCK.,0 

ISPL..FREQ. 

) 

234 

58 

4 

1     5 

1    20 

0 

1 

3 

2 

1     2 

0 

386.4 

0. 

1 

120. 

0. 

120. 

1 

0 

2 

L  . 

0. 

►       L  . 

2 

-120. 

0. 

120. 

1 

0 

2 

- 1. 

0. 

3 

-120 

V 

0.0 

-120. 

1 

0 

2 

— i . 

0. 

>     — 1  • 

4 

120 

t 

0. 

-120. 

1 

0 

2 

0. 

—  x  . 

5 

111. 

120. 

111. 

3 

4 

5 

6 

-111. 

120. 

111. 

3 

4 

5 

—  X . 

7 

-Ill 

» 

120. 

-111. 

3 

4 

5 

—i. 

—  X . 

8 

111 

1 

120. 

-111. 

3 

4 

5 

'      —  1  • 

9 

102. 

240. 

102. 

6 

7 

8 

10 

-102. 

240. 

102. 

6 

7 

8 

—  X  • 

11 

-102 

• 

2^0. 

-102. 

6 

7 

8 

>     — 1  • 

12 

102 

ft 

240. 

-102. 

6 

7 

8 

—  1  . 

13 

93. 

360. 

93. 

9 

10 

11 

1* 

-93. 

360. 

93. 

9 

10 

11 

—  X « 

15 

-93. 

360. 

-93. 

9 

10 

11 

1     —  1 . 

16 

93. 

360. 

-93. 

9 

10 

11 

I     —  1  . 

17 

84. 

480. 

84. 

12 

13 

14 

18 

-84. 

480. 

84. 

12 

13 

14 

—  X. 

•  19 

-84. 

480. 

-84. 

12 

13 

14 

20 

84. 

48  0. 

-84. 

12 

13 

14 

21 

75. 

600. 

75. 

15 

lb 

17 

22 

-75. 

600. 

75. 

15 

16 

17 

—  X « 

23 

-75. 

600. 

-75. 

15 

lo 

17 

>      —  1  . 

2* 

75. 

600. 

-75. 

15 

16 

17 

25 

66. 

720. 

66. 

18 

19 

20 

26 

-66. 

720. 

66. 

18 

19 

20 

-X . 

27 

-66. 

720. 

-66. 

18 

19 

20 

28 

66. 

720. 

-66. 

18 

19 

20 

29 

57. 

840. 

57. 

21 

22 

23 

30 

-57. 

840. 

57. 

21 

11 

23 

—  X , 

31 

-57. 

340. 

-57. 

21 

11 

23 

32 

57. 

840. 

-57. 

21 

11 

23 

33 

48. 

960. 

48. 

24 

25 

26 

3* 

-48. 

9S0. 

48. 

24 

25 

26 

—  X . 

35 

-48. 

960. 

-48. 

24 

25 

26 

36 

48. 

960. 

-48. 

24 

25 

26 

37 

39. 

1080. 

39. 

27 

28 

29 

38 

-39. 

1080. 

39. 

27 

28 

29 

—  X . 

39 

-39. 

1080. 

-39. 

27 

28 

29 

—  1 . 

40 

39. 

1G30. 

-39. 

27 

28 

29 

41 

30. 

1200. 

30. 

39 

0 

40 

0 

42 

-30. 

1200. 

30. 

39 

0 

40 

— X . 

0 

43 

-30. 

1200. 

-30. 

39 

0 

40 

-1. 

0 

44 

30. 

1200. 

-30. 

39 

0 

40 

1. 

0 

45 

30. 

1248. 

30. 

41 

30 

42 

1. 

46 

-30. 

1248. 

30. 

41 

30 

42 

-1. 

47 

-30. 

1248. 

-30. 

41 

30 

42 

™"  L  • 

48 

30. 

1248. 

-30. 

41 

30 

42 

X  • 

49 

90. 

1248. 

0. 

0 

31 

32 

1. 

50 

-90. 

1243. 

0. 

0 

31 

32 

"*  1  • 

51 

60. 

1224. 

15. 

33 

34 

35 

1. 

52 

-60. 

1224. 

15. 

33 

34 

35 

-1. 

53 

60. 

1248. 

15. 

36 

37 

38 

1. 

54 

-60. 

1248. 

15. 

36 

37 

38 

-1. 

55 

-60. 

1224. 

-15. 

33 

34 

35 

*  A  • 

»      m *  L  . 

56 

60. 

1224. 

,  -15. 

33 

34 

35 

X  « 

,     — 1 . 

57 

-60. 

1248. 

;  -15. 

36 

37 

38 

**  L  • 

»        — i  . 

58 

60. 

1248. 

•  -15. 

36 

37 

38 

J»  • 

»      —A  . 

30000000. 

.1 

-1500C 

20030. 

4. 

1 

1 

5 

1     1 

10. 

2 

2 

6 

1     1 

10. 

3 

3 

7 

1     1 

10. 

4 

4 

8 

1     1 

10. 

5 

5 

9 

1     2 

10. 

6 

6 

10 

1     2 

10. 

. 

7 

7 

11 

1     2 

10. 

„ 

8 

8 

12 

1     2 

10. 

9 

9 

13 

1     3 

10. 

^__ 

10 

10 

14 

1     3 
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APPENDIX  B 
PROGRAM  ORGANIZATION 

A.   DESCRIPTION 

The  program  organization  is  laid  out  in  the  following 
flow  charts .   The  main  driver  program  is  arranged  to  call  a 
subdriver,  XMSADT ,  and  the  optimizer  of  the  user's  choice. 
All  changes  for  replacing  an  optimizer  occur  in  MSADT.   This 
allows  for  easy  testing  of  several  optimizers  on  the  same 
problem. 

XMSADT  may  be  called  from  the  main  for  input,  analysis, 
and  output.   Printed  output  may  vary  as  the  user  requires. 
A  complete  listing  of  all  subroutines  and  their  functions  is 
given  in  Table 
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Table  XIX 
ANALYSIS  CONTROL 


ANALYSIS  ICALC=2 


ICALC=2 


NO 


CALCULATE 

OBJECTIVE, 

CONSTRAINTS 


RETURN 


YES 


CALCULATE 

SPECIFIC 

GRADIENTS 
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Table  XX 
OUTPUT  CONTROL 


OUTPUT  FOR  INFO=0/l/2 


ICALC=3 


INFO=0 


INFO=l 


(PRINT) 
-OBJECTIVE 
-G-VECTOR 


INFO=2 


(PRINT) 
-OBJECTIVE 
-G-VECTOR 
-FINITE  ELEMENT 

ANALYSIS 


(PRINT) 
-GRADIENTS 
SPECIFIED  BY 
INPUT  CONTROL 
(IDVCLC) 
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MSADT 


XMSABT 


OPTIMIZER 


Table  XXI 
PROGRAM  FLOW  CHART 

_JUNPUT) 

—  SADT4  6 


SADTQ1 


SADT0  2 


SADT11 


SADT04- 
SADT0  7 
SADT0  8 


SADT0  3 
SADT0  5 


-SADT0  6 

SADT61 SADT4  8 

-SADT0  5 


SADT1 


SADT5  0 
SADTU9 


SADT1-5 —  SADT16 
■SADT2  3 


__  SADT0  8 
— EIGRQ1 

EIGR02   EIGR03 


CANALYSIS) 


SADT7  3" 

SADT7  5 
SADT7S 
SADT7  7 


SADT36 
SADT2  0" 

SADT21 


SADT0  3 
SADT4  8 

SADT3  7 


(GRADIENTS) 


-P- 


SADT74 
SADT5  3 
SADTQ9 

SADT7  9- 
SADT17 
SADT6  2 
SADT71 


(.OUTPUT) 
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TABLE   XXII 
SUBROUTINE    USES 

FILENAME  TITLE 


EIGR01  SOLVES    EIGENVALUE    PROBLEM     | A-ALAMBD*B |     X=0 

EIGR02  SOLVES    EIGENVALUE    PROBLEM 

EIGR03  SOLVES    EIGENVALUE    PROBLEM 

SADT0  1  THIS    ROUTINE     READS    AND    PRINTS    INPUT    DATA 

AND    ORGANIZES    PSEUDO-DYNAMIC    STORAGE 
ALLOCATION 

SADT02  BUILDS    VECTORS    JC    AND    IIK    FOR    FINITE    ELEMENT 

STRUCTURAL    ANALYSIS 

SADT03  BUILDS    THE    ELEMENT    STIFFNESS    MATRIX    FOR 

TRUSS    FINITE    ELEMENT    NO.     II 

SADT04  BUILDS    ELEMENT    STIFFNESS    MATRICES    AND    ADDS 

THEM    TO    THE    GLOBAL    STIFFNESS    MATRIX 

SADT05  SUPERIMPOSES    THE    ELEMENT    STIFFNESS    VECTOR    EK 

OF    ELEMENT    1 1    ON    THE    GLOBAL    VECTOR    AK 

SADT06  BUILDS    GLOBAL    LUMPED   MASS    MATRIX 

SADT07  LU    DECOMPOSES    SYMHETRIC. POSITIVE-DEFINITE. 

SPARCE    MATRICES, THE    UPPER    TRIANGLE    OF    WHICH 
IS    STORED    IN    VECTOR    AK    WITH    LEADING    ZEROES 
OMITTED 

SADT08  FORWARD    AND    BACK    SUBSTITUTES    TO    YIELD    A    SOLN 

TO    A    SET    OF    LINEAR    EQNS ( DECOMPOSED    W/ 
SADT07) 

SADT09  PRINTS    ALL    JOINT    DISPLACEMENTS    FOR    EACH 

LOAD    CONDITION    OF    A    FINITE    ELEM.     STRUCTURE 

SADT10  PERFORMS    TRUSS    FIXED-GEOMETRY    DESIGN 

SADT11  ROUTINE    TO    ORGANIZE    ANALYSIS 

SADT14  CALCULATES       PARTIAL    DERIVATIVE    OF    TRUSS    ELEM 

K    MATRIX    WRT    COORD    DOF    11     AT    NODE    J1 

SADT15  CALCULATESS     VALUES    OF    ALL    DESIGN    AND    BEHAV- 

IORIAL    CONSTRAINTS    AS    DEFINED    BY    PROGRAM 
"SAD" 

SADT16  CALCULATES    STRESS    IN    BAR    ELEMENT    II    UNDER 

LOAD    CONDITION    JJ 

SADT17  PRINTS    FORCES    OR    STRESSES    IN    BAR    ELEMENTS 

SADT18  CALCULATES    GRADIENT    INFO    IN    FINITE    ELEM 

STRUCTURAL    ANALYSIS 

SADT19  ADDS    ELEMEMENT    MASS    AA    OF    ELEMEMENT    II 

TO    GLOBAL    MASS    MATRIX 
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SADT20       CALCULATES  THE  PARTIAL  DERIVATIVES  OF  THE 
DISPLACEMENTS  WRT  INDEPENDENT  DESIGN 
VARIABLE  II 

SADT21       CALCULATES  STRESS  AND  GRADIENT  OF  STRESS 
IN  BAR  ELEM  II  UNDER  LOAD  COND  JJ  GIVEN 
DISPLACEMENT  VECTOR  PL  AND  GRADIENTS  OF 
THE  DISPLACEMENTS  IN  VECTOR  DU 

SADT22       PERFORMS  STRESS-RATIO  DESIGN  OF  A  FIXED 
GEOMETRY  TRUSS 

SADT23       CALCULATES  WEIGHT  OF  TRUSS  GIVEN  AREAS  AND 
LENGTHS 

SADT26       CALCULATES  THE  GRADIENT  OF  WEIGHT  WRT 
RECIPROCAL  DESIGN  VARIABLES 

SADT31       CALCULATES  GRADIENT  INFORMATION  AND 
STORES  IN  ARRAY  TAY 

SADT3U       CONMIN  EXTERNAL  FOR  TRUSS  FIXED  GEOMETRY 
DESIGN 

SADT35       CALCULATES  GRADIENT  INFORMATION  AND 
STORES  IN  ARRAY  TAY 

SADT36       CALCULATES  ( XEIG-T* AM* XEIG)  FOR 

GRADIENT  CALCULATIONS  IN  FREQUENCY 
CONSTRAINTS 

SADT37       CALCULATES  EIGENVALUE  GRADIENT  INFORMATION 
IN  FINITE  ELEMENT  STRUCTURAL  ANALYSIS  AND 
DESIGN 

SADT40       CALCULATES  GRADIENTS  OF  DISPLACEMENTS 
WRT  COORDINATE  DESIGN  VARIABLE  II 

SADT41       CALCULATES  GRADIENT  INFORMATION  WRT 
COORDINATE  VARIABLES 

SADT42       CALCULATES  GRADIENT  OF  WEIGHT  WRT 
COORDINATE  DESIGN  VARIABLES 

SADT43       CALCULATES  GRADIENTS  OF  OBJECTIVE  AND 
ACTIVE  CONSTRAINTS  IN  OPTIMUM  GEOMETRY 
DIRECTION  FINDING  PROBLEM 

SADT44       PERFORMS  TRUSS  FIXED  GEOMETRY  DESIGN  AND 
EVALUATES  CONSTRAINTS 

SADT45       PERFORMS  TRUSS  OPTIMUM  GEOMETRY  DESIGN 

SADT46       READS  INPUT  INFORMATION  FOR  BAR  ELEMENTS 

SADT48       BUILDS  THE  ELEMENT  MASS  MATRIX  FOR  A 
TRUSS  FINITE  ELEMENT 

SADT49       SOLVES  REAL  EIGENVALUE  PROBLEMS  USING 
THE  SUBSPACE  ITERATION  METHOD 

SADT50       BUILDS  INITIAL  SET  OF  BASIS  VECTORS 
FOR  EIGENSOLUTION  BY  REDUCED  BASIS 
METHOD 

SADT51       CALCULATES  STRESS  AND  GRADIENT  OF  STRESS 
IN  BAR  ELEMENT  II  UNDER  LOAD  CONDITION 
JJ,  GIVEN  DISLPLACEMENTS  CONTAINED  IN 
VECTOR  PL  AND  GRADIENT  OF  DISPLACEMENTS 
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SADT52 

SADT53 
SADT54 

SADT6  1 
SADT62 

SADT71 
SADT73 
SADT74 
SADT75 

SADT76 

SADT77 
SADT79 
MSADT 

XMSADT 


IN    VECTOR    DJ 

PERFORMS    1-DIMENSIONAL    SEARCH    IN 
COORDINATE    DESIGN    SPACE 

PRINTS    MEMBER    INFORMATION    FOR    BAR    ELEM 

CALCULATES    AND    PRINTS    CENTER-OF-GR A VITY 
AND    INERTIAL    PROPERTIES    OF    FINITE    ELEM 
STRUCTURE 

BUILDS    GLOBAL    MASS    MATRIX 

PRINTS    NEIG     EIGENVALUES    STORED    IN 
EIGVAL,    AND    THEIR    CORRESPONDING 
EIGENVECTORS     STORED    IN    XEIG 

PRINTS    G    VECTOR    OF    CONSTRAINTS 

CALCULATES    GRADIENT    INFORMATION 

PRINTS    TAY    ARRAY    OF    GRADIENTS 

CALCULATES    GRADIENT    OF    WEIGHT    WITH 
RESPECT    TO    AREA    DESIGN    VARIABLES 

CALCULATES    GRADIENT    OF    WEIGHT    WITH 
RESPECT    TO    COORDINATE    DESIGN    VARIABLES 

CALCULATES    GRADIENT    INFORMATION 

PRINTS    COORDINATE    INFORMATION 

DRIVER    PROGRAM    FOR    USING    THE    ABOVE 
SUBROUTINES.     MSADT    MAY    BE    COUPLED 
TO    OPTIMIZER    OF    USER'S    CHOICE. 

SUBDRIVER    PROGRAM    FOR    COUPLING 
SADT    ROUTINES    TO    MSADT 
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